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Abstract 

New results for the fermion propagator in a laser background are presented. We show that the 
all orders electron propagator can be written in a compact and appealing form as a sum of sideband 
poles with a matrix wave function renormalisation and a matrix valued mass shift. This last result 
is essential in the fermionic theory if we are to maintain that both the mass and its square pick up 
a correction only at order e^. A perturbative verification of our results is carried out. 

PACS numbers: 11.15.Bt,12.20.Ds,13.40.Dk 


1 


I. INTRODUCTION 


In this paper we will discuss fermionic propagation in background fields. Specifically we 


are interested in electrons propagating in a plane wave backCTound |l|. This is important for 
understanding recent work with high power laser facilities |2i|3|. The result we will obtain 
is simple and has a clear physical interpretation. A surprise is that all the renormalisation 
factors are matrices. 

The propagator is one of the basic building blocks linking physical processes with the 
structure of quantum held theory. We recall that for a scalar held, 4>{x), the Feynman 
propagator is given by the time ordered product 


{0|T.#.(i)#(j/)|0) = , (1) 

J + te 

where the helds are taken to be free and the bar in the measure indicates appropriate 
factors of 2n. In the presence of loop corrections the form of the propagator is preserved 
and typically becomes 


{O|T0(a:)0t(l/)|O) = d^p - 


iZo 


( 2 ) 


p^ — {m? + 5m?) + ie ’ 

which can be interpreted as the propagation of renormalised helds with a renormalised mass. 

A similar, but richer, structure in the propagator is found when the scalar held interacts 
with a laser plane wave background jd] jsljb] jrljs] j9|[w|[ll|[^| 13| ^ 


OO « 


-ip-{x-y)_ 


iZ. 


(n) 


{p + nk)‘^ — (m^ + 5m^) + ie ’ 


( 3 ) 


where we see a sum over so-called sideband states 


I^[l6|[l7| and note that the mass shift 


is common to all sidebands. The exact form of the mass and the momentum dependent 
wave function renormalisation factors depend on the details of the laser background. For a 
linearly polarised background, the propagator has the form 


E 


f.4 -rHx-y) 

’ ^ {p + nk)^ -ml + ie' 


( 4 ) 


where Up and Vp are laser parameters that we dehne in (1201) . The mass shift is here given 
by ml = m? — \e^a? where is the (spacelike) laser amplitude and this mass shift does 
not acquire higher order corrections. The wave function renormalisation for the various 
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sidebands is written in terms of a prodnct of generalised Bessel fnnctions which can be 
expressed in terms of normal Bessel fnnctions as 

OO 

Jfi(e'Up, e Vpj ^ ^ ’Jn—2r(^(^Up) ^p) • (^) 

r=—oo 


The hne details in (jl]) are polarisation dependent bnt the broad strnctnres for a plane wave 
laser backgronnd, inclnding the snm over sidebands and the mass shift being at order 
only, are common. 

Note that in the literatnre, the two point fnnction in a laser backgronnd is widely referred 
to as the propagator. However, throngh interactions with the backgronnd, this two point 
fnnction inclndes diagrams where the initial and hnal momenta of the matter held are 
distinct. In this paper, and in onr snmmary above, the focns is on the diagonal part of the 
two point fnnction which does not inclnde momentnm altering vertex effects. This is what 
we call the propagator in this paper and this is given for the scalar theory in (HI) . 

As in onr previons work, we are gnided below by experience from infra-red physics. A 
charge propagating in a laser is indistingnishable from a charge which absorbs and emits 
the same nnmber of laser photons. These laser indnced degeneracies parallel the soft and 
collinear degeneracies of the infra-red regime in QED and QCD 18| 19| ^ 211. Under¬ 
standing the mass shift in a laser backgronnd may also help to clarify the cnrrent versns 
constitnent mass distinction in QCD 221. 

Tnrning now to fermionic QED, the pertnrbative propagator withont a backgronnd has 
the form 

'■ iZ2 

nl 4- ze ’ 

( 6 ) 


(0|TV>(a;)V>(z/)|0) = 

= [d'^pe 


— {m + 6m) + ie ’ 

-ip-ix-y) iZ2{f + m + 6m) 

p^ — {rn + 6m)‘^ -|- ie ’ 


where we note the presence of single powers of m. 


The extension o 
sidered jd, 5, 


this theory to plane wave laser backgronnds has been previonsly con- 


231, [2^. The resnlts obtained also display the sideband strnctnre of the scalar 
theory bnt the two point fnnction is qnite involved, see, e.g., (5.7) and (5.8) in Ref. jd], and 
does not clearly resemble a renormalised propagator. In particnlar in the nnmerator of the 
propagator obtained by Reiss and Eberly there are two gamma matrix strnctnres, one 
proportional to ^ -f m and one proportional to the laser momentnm, each mnlti- 
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plied by combinations of Bessel functions. Similarly, additional gamma matrix structures 
are obtained in Appendix A of Brown and Kibble js] and by Ritus jn |. 

The mass shift in the fermionic theory also satishes the spin independent relation % [S] 


2 i ir\2 2 122 

= [m + dm) = m — -^e a , 


(7) 


to all orders in the coupling. This might naturally suggest, see for example section 40 of 251 


and more recently 


261, that 


m* = m + 6m = m\ 1 


e^a? 
2m^ ’ 


( 8 ) 


which would, however, mean that 6m would acquire corrections to all orders in the coupling 
through an expansion of the above. 

Below we will first show that an explicit perturbative calculation disagrees with this 
expectation for the mass shift. We will then present an all orders calculation of the fermion 
propagator which will reveal a simple form which can be interpreted, as for the scalar theory, 
in terms of mass and wave function renormalisation factors. An unexpected feature is that 
all the renormalisation factors we hnd, including the mass shift, are matrices. The mass shift 
will itself only be corrected at order e^, but still satisfy ([7]) to all orders. As this is surprising 
we will independently verify all of the structures we obtain via a low order perturbative 
calculation. Finally we will present some conclusions. Some technical details are presented 
in appendices. 


II. PERTURBATIVE EXPANSION 

The Lagrangian density describing the interactions of a fermion with the laser is 

C = — m)'4). (9) 

In this the covariant derivative is given by + ieA^ and, in a linearly polarised 

background, 

Af,{x) = af,cos{k ■ x), ( 10 ) 

where k = fco(l, 0, 0,1) is the null vector characterising the laser. Note that the gauge fixing 
condition used is /c • a = 0. 
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FIG. 1. The three point Feynman rules for linear polarisation 
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FIG. 2. Leading order contributions to the propagator 


The Feynman rnles for the vertices in this linearly polarised theory are given in Fig. 1. 
For this polarisation, the rnles are identical for the absorption and the emission of a laser 
photon. The two diagrams that contribnte at lowest order in the conpling are shown in 
Fig. 2. The hrst diagram yields after a partial fraction expansion 

^ + m) + (p + kf - 


and the second diagram is obtained by replacing k —)■ —k. In our perturbative calculations 
we suppress the ie prescription. The gamma matrix structures are now expanded depending 
upon whether the denominator contains poles in or is in an adjacent sideband. In 

this initial discussion, we shall restrict ourselves to the central sideband alone. In the double 
pole coefficient the 1^ term alone survives when the two diagrams are added, while in the 
single pole coefficient the ^ term cancels. Thus the two diagrams contribute to the central 
sideband 

^ 1 ^p-k 


le 


+ m 


8{p-ky 

In the first term we rewrite 


+ m 


+ m) 


p2 _ ^2 


+ [f + m 


+ m) 


(p2 _ 17^2^2 


( 12 ) 


{f + m)(j{f + + m) = {f + m)^ + ‘2^p-a{f + m)) . (13) 

The hrst part here will not contribute any pole to the central sideband and we have verihed 
that the sum of all such terms without poles cancel when the contributions from all sidebands 
are summed. We will therefore just drop them below. The double pole term in flT^ will also 
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generate single pole terms in a similar fashion. The sum of all the terms containing poles in 
the central sideband is found to be 


f {p-aY f + m 1 I ^ f + ^ \ 

\ 2{p-kY p"^ — m? Ap-k p"^ — rn? 2^ (jp'— w?Y ) 


(14) 


We should now compare this with the general form of (El). At this leading order we expect 
three contributions. Writing = 1 + (where SZ^^^ is expected to have a matrix 

form) this would give 


le 


SzY\p + m) 

p2 _ ^2 


+ 6m— 
P- 



1 2 1 \ 

2 (p2 — m^)2y ’ 


(15) 


where we have used d?]) in the last term. 

In 0141) the double pole term is the expected mass shift from the expansion of the denom¬ 
inator of the central sideband. The hrst term gives us information about the wave function 
renormalisation factor. However, the second term in (ITT)) , which involves aP as would be 
expected of a mass shift, contains a |^-factor and is not obtained by simply expanding the 
numerator in ([6]). We will therefore now carry out an all orders calculation of the propagator 
to resolve this puzzle and also determine the form of z'YK 


III. VOLKOV FERMION 


We now want to construct the spinor held which satishes 


{ilp — m)'ipy{x) = 0 . 


(16) 


This Volkov equation (1161) is solved by 


Yyix) = Id^P^i J^i{x,p)DYx,p)Ul°‘\p)a^Y^{p) 


+ Teix, -p)DYx, -pWYYpYI^'^Xp) 


(17) 


There is a lot of structure to this solution that we need to unpack carefully. 

The on-shell momentum = {E*,p) satishes the shifted on-shell condition of the quasi¬ 
momentum: 

p^-ml = E;^-f-ml = 0. (18) 
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Just as for the scalar case, plane waves are distorted in this background and become equal 


to 


D (^X p) = e-*P'^ehe“pSin(fc-a;)+e2DpSin(2fc-x)) 


= 




( 19 ) 


In this expression we have dehned the laser variables Up and Vp by 

2 


p-a a 

Up = and Vp = 

p-k Sp-k 


( 20 ) 


The rest of equation (1T7|1 reflects the fermionic nature of QED. Focusing initially on the 
hrst term in we see that in the Volkov equation the derivative will only act on the 

factors J^i{x,p)D^{x,p). The prefactor J^£{x,p) is a matrix given by 

Mix) 


Tiix.p) = 1 + e 


‘Ip ■ k 


( 21 ) 


The action of the covariant derivative is then given by 

il^J^iix, p)D^ (x, p) = Tiix, p)D^ (x, p)ip + 2e^vpp). (22) 

Similarly, in the second part of flTTl) we get 

iJpJ^iix, -p)D^ix, -p) = J=e,ix, -p)D^ix, -p)i-p - 2e^Vpp). ( 23 ) 

The Volkov spinors are then introduced to satisfy 

if + 2e^Vpp)U^'^ ip) = m ip ), ( 24 ) 

and 

if + 2e^Vpf)V)f> ip) = -m ip). ( 25 ) 

The precise form for these Volkov spinors is presented in Appendix A. 

Combining these results we see that (x) as dehned by equation (ITTI) does indeed satisfy 
the Volkov equation (ITHIl . 

Having constructed the Volkov held, '^^(x), we now simply state the expression for the 
Dirac adjoint 


= d^q^\Up[q)F,(y,-q)DHy,q)aW\q^ 


+ V‘f '’(q)Ti(y,q)Dl(y, -q)b[d\q) 


( 26 ) 
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Returning to the on-shell condition fl2^ for the Volkov spinors, we can directly reinterpret 
this as 

(p) = (m - 2e^vpf:) (p), (27) 

that is, in the spinor theory, the mass shift has a matrix form with m* = m + Srfi with 

2 2 

5ffi = —2e^Vp^ = — . , ^ ■ (28) 

4p-fc 

Let us now show that this is consistent with the result familiar from the scalar theory. We 
multiply both sides of (HT} by p and anticommute the p factor through the p, to obtain 

p2 ^ _ 4g2^_p . ^ _|_ 2e^Vppp) (p). (29) 

We now reuse fl27)l in the hnal term and, recalling that = 0, obtain the spin-independent 
on-shell result 

p^ (p) = ml (p) = (m^ - \e^a^) (p). (30) 

Note that this requires 6ffi to be a matrix, as one would otherwise generate corrections to 
at order e^. 

However, we cannot yet identify 6m in (1T5|1 from this, as if we were to replace 6m by 6ffi 
we would obtain (1141) but with the wrong sign for the $ term. We shall though see later 
that the result for 6ffi is correct and that the naive expectation (IT^ is wrong. 

IV. THE PROPAGATOR 

We are now in a position to calculate the Volkov propagator 
v(0|T?/>y(x)V!y(p)|0)^ = e{x^ - y^) ^{0\i/j^{x)^^{y)\0)^ 

- 0(p° -x°)y(0|V;y(p)V’v(a^)|0)v , 

where the Volkov vacuum is identihed through the modes in fll7p and fl2^ . Note that the 
creation and annihilation operators satisfy 

{«1“’(P). “fb)} = (27r)=Ri“'’i'='>(p - 9). (32) 

Looking at the hrst vacuum expectation value we have 

v(0|^v(^)^v(2/)|0) = ^J^eix,P)ip + m - 6Tfi)J^e{y, -p)D^{x,p)Dl{y,p). (33) 




In a similar way we obtain 


v(0|V'v(2/)^v(a^)|0)v = -p){p-m - 6rfi)Te{y,p)D^{x, -p)D\{y, -p). (34) 

We now use these to build up the propagator. To this end we require the following results 
which follow from flT^ 

T,(x,p)D,(x,p) = ^ (35) 

n=—oo ^ 

where we have condensed our notation for the generalised Bessel function so that 

Jn(£''Up, e Up) , (36) 

and 

J^{p) = J^_,{p) + Jx,+i{p). (37) 

In a similar way we obtain 

T,{x,-p)D,(x,-p)^ f; e*V"‘-’(j„(-!5) + e^J„'"(-p)), (38) 

n=—oo ^ 

while 

Ti(v,-p)Dl(y,p)= f; e'”e-"‘'»(j-„(p) + e^J^(p)). (39) 

n=—oo ^ 

and 

Hv,p)Dl(y,-p)= f; e-»»e-’“-''(j;.(-p) + e^J?(-p)). (40) 

n=—oo ^ 

Inserting these expansions into fl33l) and fIMl) and using the decomposition 


OO 

^v(^) = '^n{x) , 

n=—oo 


(41) 


gives the propagator dehned through the diagonal contributions to the two point function 

1 


A0\Mx)My)\0), = J„(p) + 


2E* 


ip-k 


X ip + m - Srfi)(^J„{p) . 


(42) 


and 




2E. 


Ap-k 


X{f-m-5Tfl) [jn{-p) + ^^^'^ni-p)) ■ 


( 43 ) 
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The Volkov propagator then becomes equal to the sum over n of 


/ 


'dp'^ d^p i 


+nk°){x°-y °) i(p-nk)■ [x-y) 


271 2E* — ie 



( 44 ) 


^-iip°-E*-nk°)(x°-y°)^-i{p+n!£i-{x-y) 


X (Jn(-P) + ip-m- 6lffl) {Jni-P) + 6 


— m — 



Combining the variables and p® into an off shell four moment p^ (details of which are 
given in Appendix B) we can write the fermionic Volkov propagator as 



X 


i{f + np + m — 5rpi) 


(45) 


(p + nky — ml + ie 


X (Jn(p + nk) + ^^^J'ni.P + ■ 


Note that here we have removed the bar over the mass shift matrix dehned earlier in fl28|l 
since the momentum in (1451) is now off shell. 

V. INTERPRETATION 

The naive expectation might have been that we would hnd a result of the form 



(46) 


however, what we have found in fl45p looks very different. It contains matrix dependent 
wave function renormalisation constants on both sides of the usual gamma matrix structure 
expected in a fermionic propagator as well as a minus sign in the mass shift, 5rpt in the 
numerator of (I45|) which we will now explain. 

Setting u = 0 for simplicity and ignoring the ie, we make the simple observation that 
because the mass shift is a matrix, then it can naturally be regrouped as 


1 


1 


(47) 


f — {m + 5rpi) — difi) 


— m 
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which means when we carry out the standard algebraic rearrangement we obtain 


1 f — 5^ + m 

— Sifi) — m {f — 8ifiY ~ ' 

f + m — brja 
jp' + e^a^/2 — ra^ ’ 
^ + m — bija 

p2 _ ^2 


( 48 ) 


Here we recognise the standard mass shift in the denominator and the minus sign present 
in the numerator of 

We finally define the wave function renormalisation matrix by 



Jn{p + nk) + + nk ), 

Ap-k 


(49) 


and note that 


:= V = Jnip + nk) + e^JniP + nk). (50) 

This yields a very appealing way to write the sideband propagator. If we reinstate the sum 
over sidebands, we see that the Volkov propagator for fermions can be written as 

OO OO « / . - - 


^{0\Tpn{x)pn{y)\0)^ 

n=—oo 




n=—oo ' 


d pe 


f + nlf: — {m + 5ifi) + ie 


r(n) 


(51) 


VI. DOES THIS AGREE WITH PERTURBATION THEORY? 

We will now check that our result for the propagator fl51l) agrees with perturbation 
theory by calculating the lowest order corrections to the free propagator extending our 
earlier result fllTll to the adjacent sidebands. 

We will require the expansions of the generalised Bessel functions 

(PvP 

Jo{p) = Jo(eUp, e^Vp) = 1-+ 0{e ^), 

4 (52) 

PH 

J±i{p) = J±i(eMp, e^Vp) = + 0{e^). 

We should immediately note that fl5T]) is indeed an expansion in e^. In fl45l) we see that for any 
value of n the bracketed terms containing the combinations of generalised Bessel functions 
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are either both odd or both even in e due to the property J„(—x,?/) = (—l)"'Jn(x,?/). As 
there are two bracketed terms, the overall propagator will be an expansion in e^. 

Let us hrst consider the central sideband at leading order. The wave function renormal¬ 
isation factors are 

= \/if = (53) 

as the matrix factors cancel at this order here. Using it is now simple to reproduce all 
of flIT|) from an expansion of (l5Tl) . This calculation is sensitive to both the matrix form of 
the mass shift and the sign in the numerator of fH8|) . 

In the hrst upper sideband we proceed similarly to flT^ but adding and subtracting 
^ terms as follows 


{f -h m)(fi{f + lli + + m) = {f + + Ij: + ^ m - f), (54) 


which may also be repeated on the left. This generates zero pole terms that cancel with 
those found in the central sideband and from the hrst lower sideband. The terms containing 
poles in the hrst upper sideband are thus found to be 


le 



f + $ + m 
m? 
2 


{p + ky 


1 p-a 

+ 


p-a 

p-k 


Ap-k^ {p + ky 
1 




{p ky 






(55) 




and similarly with k ^ —k for the hrst lower sideband. (The adjacent sidebands will hrst 
acquire mass shifts at the next order in the coupling.) 

This calculation also agrees with flHTll at this order. In the hrst upper sideband, we note 


that 





(56) 


at leading order in e. There is a similar expression for y . Substituting this into fl5T|) 
yields fl55l) . The hnal a^-dependent term in fl55l) may appear surprising in this sideband as 
there is no contribution from the mass shift in the perturbative expansion at this order. It 
arises via a product of two ^-factors in the gamma matrix structures. We note that this 
calculation is sensitive to the matrix structure of the wave function renormalisation factors. 
The hrst lower sideband is again simply reproduced by changing k —)■ —k everywhere in fl5^ . 
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As the agreement across these three sidebands depends upon the form of the various 
structures in flSTl) . it verihes the emergence of a matrix mass shift and matrix wave function 
structures ordered to the left and right. 

We also comment here that the propagator in the laser background cannot be written in 
the form 


(57) 


— m — S ’ 

with —iS being the self-energy. This is because the usual derivation of fl57l) involves the 
identihcation of the perturbative expansion 


— m 


+ 


— m 


-(-ffi)- 


— m 


+ ••• 


(58) 


as a Taylor expansion of (j57|l . However, this requires that f and S commute. This is 
true in normal perturbative QED since the self-energy then only depends on f and scalars. 
However, in the laser context, there are 1^ and (ji factors which do not commute with ^ and 
hence obstruct an expression of the form fl57|l . 


VII. CONCLUSIONS 


In this paper we have derived an expression for the fermion propagator in a linearly 
polarised plane wave. In QED, as well as in the scalar theory, all of the laser induced 
structures in the propagator can be understood as multiplicative renormalisations of the 
free propagator. The distinctive feature highlighted in this paper is that for fermionic QED 
these multiplicative factors are matrices. Matrix structures have already been seen in earlier 


work (^[5|][lll], but here we have seen that the electron mass shift is a matrix such that 
both the mass and also the mass squared only pick up corrections at leading order in the 
coupling. Although we have focused in this paper on linearly polarised lasers, our approach 


can be generalised to other polarisations, see 


13| . We expect that such a matrix valued 


mass shift will persist in more realistic models of laser pulses 27|]. We note that matrix 
wave function renormalisation constants have been seen previously in studies of the infra¬ 


red 


although the mass shift did not there involve a matrix. We have also checked our 


all orders calculations through a hrst order perturbative expansion of the propagator which 
tested all of the key features of our overall result flSTj) . 
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The results presented in this paper will be important for calculating scattering processes in 


where the propagator is used to ca 


QED in laser backgrounds, see for example 
ble Compton scattering, and also 30[|Mi for this process. More generally see 321 


culate dou- 


^| 35 | 


for other processes. It is intriguing to speculate how the matrix renormalisation constants 
impact upon the calculation of S-matrix elements and the physical predictions of QED in a 
laser background. 


Appendix A: Volkov Spinors 


The defining equations fl2TD and allow us to construct the Volkov spinors by suitably 
boosting the static spinors 




fl\ 
0 
0 

W 

To this end, we define 




/o\ 

1 

0 

VO/ 


, v«(0) = 


/o\ 

0 

1 

VO/ 


, v^^Vo) = 


/o\ 

0 

0 

VV 


^2m{m + E* + 2e‘^Vpko) 


and 


= 


'2m{m + E* + 2e‘^Vpko) 
With these conventions we have the inner products 


= and (P) = 


(Al) 


+ 2e\^ + m)W(A(0), (A2) 


:(-^-2e%|^ + m)VA)(0). (A3) 


(A4) 


as well as the tensor products 




P + 2e^Vpp, + m 


2m 


and 


vpiPW'ip )= 


(AS) 


(A6) 
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Appendix B: Volkov propagator details 


The argument taking us from the initial expression for the Volkov propagator in Equa¬ 
tion (HT|) to the final covariant expression fHSD involves several key steps that we outline 
below. 

Starting from fHT|) we make the substitution ^ E* — — nk^ in the hrst part and 

p^ ^ p^ + E* + nk^ in the second to give 



f _!_ Q-ip° -y°) gi(p-nk)-{x-y) 

\E* — p^ — nk^ — ie 

X (Jn(p) + ip + m- 5rfl) {Jnip) + 

_!_ Q-ip°{^°-y°)Q-i{p+nk)-{x-y) 

E* + p^ + nkP — ie 

X {Jn{-P) + e^jt{-p)){p-m - bTfl){jn{-p) + ^^J'ni-P))) ■ 


(Bl) 


If we now shift the three momentum in the first part by p ^ p -|- nk and by p —)■ —p — nk 
in the second (noting a change in sign in 6rfi) then we have 


d^p 


o-ip-{x-y) 


2E, 


p-\-nk 


- po - nko - ie 


-{Jnip + nk) e 


4(p -|- nk)-k 


np + m — Srfi) {Jn{p + nk) + e 


4(p -|- nk)-k 


JniP + nk)) 


JniP + nk)) 


E;+nk + 


-{Jn{p + nk) + e 


4(p -|- nk)-k 


-Jn{P + nk)) 


{-f + np - m + 5rfi) {Jn{p + nk) + e 


4(p -|- nk)-k 


-JniP + nl!)) 


(B2) 


In this we have simply shifted on-shell momenta in the hrst part but in the second part 
noted that under the shift 


p -)■ {Ep+nk, -P - nk ), 

( nE ), (B3) 

= —p -|- nk , 

where we denote by the tilde the negative energy on-shell four-momenta. 
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Combining these terms over a common denominator allows us to extract the covariant 
propagator and we get 


+ + m — 5ffi) {Jn{p + nk) + e 


- (£^p+nfc - - nk^) {Jnij) + nk) + e 


4(p + nk)-k 


Jn{p + nk)) 


4(p + n/c)-/c 


-J^{p + nk)) 


{—f + n$ -m + 5rfi) {Jn{p + nk) + e 


4(p + nk)-k 


-Jn(P + nk))) . 


(B4) 


We now recognise the two terms as simply the residues associated with the two poles of the 
covariant denominator. Hence we recover the covariant expression (05]). 


Acknowledgements: We thank Tom Heinzl, Anton Ilderton and Ben King for helpful 
discussions. 
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